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the local cooling time and the thermal instability timescale. This model was subsequently 1° the present paper, we will adopt the procedure of TB and reanalyse the data in the 

studied numerically where it was found that the flow structure was basically the same for following manner. We will first assume a multiphase model with which we can compute 

the two choices of timescale and that the resulting distribution of accreted material was the surface brightnesses as seen by the HRI and IPC. We can, furthermore, compute the 

in general steeper than an isothermal mass distribution. When the cooling time model hne fluxes of the lines seen by the FPCS and the mass profile after assuming hydrostatic 

was applied to data from M87 and NGC 1275, it was found that the model provided a equilibrium. The parameters of the model are then adjusted to fij the data. We will then 
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expected to be the coolest (Canizares et al. 1982; Lea et al. 1982), using the higher, field column density without a cooling component does not result in a better fit of the models to 

of view averaged temperature may ascribe too much of the low temperature emission to the SSS data. In this paper, we will be able to consider this question in the more general 

the cooling component of the gas. context of the two previously mentioned cooling flow models. This wifi allow us to see 


whether the addition*] spectroscopic and imaging data besides that from the SSS that we w ' begin by stating our assumptions. First, we will suppose that the thermal emission 

will include in our analysis can be consistently explained along with the SSS data when bom the intracluster gas around M87 comes from two distinct components. The first 

the column density of absorbing gas is taken to be larger than the Galactic value. This component, or phase is an ambient, non-cooling distribution of hot gas that is taken 

will be done by performing fits to the data assuming two possibilities for the distribution 10 be spherically symmetric, in ionisation equilibrium, and describable by density and 
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section IV. We will present our derived mass accretion profiles in section V and we will ( 1989 '. P nTate communication), to compute the X-ray emission. We will not discuss the 

draw our conclusions in section VI. validity of the other assumptions further as this issue is considered in the work of White 

and Saraxin. In addition, the assumptions that have adopted with regard to the ambient, 
IX. METHOD OF SOLUTION non-cooling, medium have been discussed in section III of TB. 


We now further assume, as in TB, that the density and temperature profiles of the where q is & parameter that determines the " efficiency” with which the cooling time mass 

ambient medium are parameterized by the forms drop-out occurs. Here, we have the definition 
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In performing our fits to the d»t», we will foUow the treatment of TB with respect other parameters associated with the temperature profile to vary with the fit. It will turn 

to the elemental abundances. We will first fix the relative abundances of the metals other out that we will have to modify this part of our procedure below because the resulting 

than oxygen to be the solar values (Allen 1973). We will then fix the overall abundance temperature profiles from a fit with a fixed T* is inconsistent with spectral data from the 

to an initial guess and perform the fit to the data. The equivalent width of the 7 keV IPC. This aspect of the fit constitutes a departure from the procedure of TB which we will 
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our density profile as parameterized in eq. (1) will be -<*! - a 3 ) which must then equal re-evaluate the evidence for the presence of additional absorbing gas in the cluster. This 

the fixed slope as determined from the IPC surface brightness profile, which is 1.36. In will be done by first assuming that the value for the absorbing column density quoted by 

addition, since that temperature of the gas at large radius is not well determined by the White et al. (1991) applies uniformly across the duster. We will then determine if the fit 

available data (see section II of TB), we will first fix X*, to be 3.8 x 10 7 °K and allow the to the low energy part of the SSS is correspondingly improved over the case where only 
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In Figure 8, we .how the mass profile, that corre.pond to our model, when the «- to x 1 due to the FPCS line, and tho~ due to the IPC part of the data indicate a fairly 

.umption of hydrostatic equilibrium is adopted. These are plotted along with optical good fit to the data. However, the fit to the mass and HR] parts indicate a poor fit. We 

determinations of the mass by a series of authors. As we mentioned in the previous sec- now first consider the IPC and FPCS parts of the fit. In the case of the IPC, the data 

lion, we indicate the mas. determination, of Sargent et a). (1978) from spectroscopic are relatively straightforward to understand and use in the following sense. The source. 



of emission other than that due to the thermal bremsstrahlung emission are relatively of the residuals with radius is basically the same with the largest residuals arising from 

absent in this part of the cluster. Furthermore, point sources have been removed and the the second and third outermost points. The innermost point contributes the next largest 

circularity of the X-ray isophotes indicate that the assumption of spherically symmetry is amount to X 2 in that the emission is somewhat underproduced at this point. Now, the first 

probably acceptable. Furthermore, the error bars are relatively well understood and can be issue that must be considered here is the manner in which we have treated the data for this 

used directly in a X 2 fit. We conclude therefore that the value of X 7 gives good statistical instrument, which was discussed extensively in section Ill.(a) of. TB. We noted that the 
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case, it was a found that an acceptable fit to the entire X-ray data set was obtained. We Since the background subtraction is performed by subtracting a constant level from the 

can therefore compare the residuals in the HR1 part of the fit for the acceptable single phase data, this could imply some correlation among the resulting data points, especially in the 

fit (see Figure 12 of TB) and the multiphase fit. From this we see that the distribution outermost points which have values that are not much larger than the background level 


itself. The value of X 7 would then suggest a worie fit than is actually the caae. This is We now tr X understand the effect on the single phase fit of including a cooling 

also a consideration which may apply to the IPC. Note that for this instrument additional component. First, we see that the rather significant decrease in the contribution to 

errors were added to the statistical errors in order to account for background subtraction 01 the FPCS brought about by the inclusion of a cooling component arises out of a 

and vignetting uncertainties. Furthermore, the data were binned into 1* rings, but the 8 ener,J decre “ e ^ overall X 2 in the multiphase case but offset by small increases in the 
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including the emiaakm due to the power Uw component. Now, by comparing this figure to fit in all cases implies that the the density and temperature profiles will be so adjusted 

Figure 21; of TIB, whiqh shows the single phase fit to the SSS spectrum, we see that thei the in the new fit such that the total emission that comes out in the Fe L part of the SSS 

inclusion of a cooling component does not improve nor greatly cbangq the fit. As in the spectrum must be comparable from one case to the next. 



b) Results with addition*] column density The FPCS residuals are given in panel (c) of Figure 14, and the predicted mass profile is 

given as the solid lines in Figure 15. 

In this section we present our results for the fit to the data assuming the cooling time 

model but with a value for the absorbing column density larger the Galactic value. We We see immediately from these figures that the major improvement in the degree of 
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For convenience, we now first describe the contents of the next few figures and then additional fits with larger assumed values for the abundance of iron, w C find that indeed 

discuss them simultaneously. In figures 12 and 13, we have plot as the solid lines the we must again increase the abundance of iron to at least 2 times the solar value (and 

predicted HR1 and IPC surface brightness profiles, respectively. Again, the other lines oxygen to at least 24 times the solar value) in order that the predicted equivalent width 

corresponds to fits where the column density varies with radius and will be discussed later. come to within 1 cr of the measured value. Therefore, this fit seems untenable in the same 
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TB CM! be referred to for the details. Now, the increased emission at around 0.65 keV in not occurred, which we can see by consulting Figure 17. Here we plot the flux of the eight 

the SSS spectrum is due to the mcreased oxygen abundance required to fit the data with observed FPCS lines as a function of temperature in the case where we have fit the data 

the presence of additional absorption. This arrises in the fit because of the FPCS data on assuming Galactic absorption (the nominal case). Note that this plot is the cooling time 
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standpoint, the FPCS and the SSS have different apertures, but they are both strongly mass profile). We have furthermore fixed the abundances of all metals to be half solar and 

weighted toward, emission in the central regions of the cluster. A situation where much the column density of absorbing gas to be the Galactic value. We have done this with the 

of the emission from the cool lines as seen by the FPCS onginates outside of the SSS field following in mind. The large amount of emission in 0V1II Lyo seen by the FPCS is the 

of view would imply that the cool gas resides outside of about 3‘ from the center of the dement which drives up the abundance of oxygen relative to the other metals and which 


also drives up the emission at 0.65 keV in the predicted SSS spectrum. By holding the SSS spectrum then goes far above the data. That this happens is easy to understand. The 

abundance of oxygen down to the same level as the other elements and by removing the significantly decreased emission at low energies again requires the emission at high energy 

agent that pushes for the higher oxygen abundance, we can assess whether we would still to he enhanced so that the HRI data can be well fit. The reason that the emission at higher 

hare as significant & problem with the overprediction at 0.65 keV as before. In Figure 18 energies is so much overproduced in this case is that the HR! effective area is higher at the 

we show the predicted SSS spectrum for this fit (light solid histogram) along with the data. lower energies of the OVIII Lyo line. Therefore, to make up for Jthe loat emission at low 
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wumed that the abundance of all the elements is half the solar value, but this time we from various sources. Spectral work from the IPC (Fabricant et al 1980; Fabricant and 

have: assumed that the total column density is 1.75 x 10* W. We see that we now have Gorenstein 1983} indicate that the column density in the outer regions away from M87 

’W good agreement at the lowest energies and up to about 1.1 keV where the predicted is consistent with the Galactic value. On the other hand, work on data from the SSS 


(Lea et al. 1982, White et aJ. 1991) have always indicated a higher value for the cohimn particular values will be made clear below. We note that as with our other fits, we have 

density. Since the SSS observations were taken of the central regions of M87, the two again specified that the abundance of metals other than oxygen is fixed at 0.7 times solar 

disparate values could be consistent if the column density were taken to be radially varying and we have allowed the abundance of oxygen to vary. We now first discuss the quality of 

with a higher value in the center of the cluster and the lower Galactic value further out the fits to the data. In Figures 10 and 11 we show the density and temperature profiles 
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u for previous cues. Here, panel (d) give, the cue where t n „ i. 10'« year., panel (e) *• much *"»dler than the Galactic value. Only at the innermoat region, if it comparable 

correspond, to 5 X 10* years, and panel (f) correspond, to 10’ year.. Finally, the maaa to the G »lactic vlJu '' We c “ now undectand the effect of this extra absorption on the 

profiles are shown in Figure 15. The predicted profile, are again .o similar that they lie on modcl Predictions for each instrument. First, .ince the IPC look, at the part of the data 

top of each other along the light solid line. By comparing the predictions in the cue. with ° ut * ide ° f * bout 3’, it will not be effected. A. for the FPCS, about half of its emission in 

radially varying column density with the prediction, of the cue with Galactic absorption, »*>' various “me. out at radii outside of 3' and the other half-come, from inride of this 
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10® years- We note first the reasonable property that the shorter the assumed lifetime of plausible from the standpoint that the maximum radius lies within the bounds of the part 

the cluster, the smaller the additional extra column density, We can now understand why of the cluster that we have modelled. Since our modelling has only been done for gas that 

the fit is not changed by including extra absorption. For most radii, the column density w relatively close to the center of M87, if is rather large, and if the contribution to 
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a fit with a flattened power law to assess possible improvements in the fit; see below). 

We finally turn to the predicted SSS spectrum. We have not shown our spectrum here 

since it is again basically the same as that pictured in Figure 9. The previously considered In Figures 20 and 21, we plot as the solid line the density and temperature profiles 

problems with the overprediction at 0.65 keV and underprediction at around 1 keV persist that result from our fit to the data assuming a rou j t = 3.5. The dotted lines in these figures 

even in our present fits. This is understandable from the standpoint of our discussion and the subsequent two figures correspond to a fit assuming extra absorption which we will 
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which is comparable to the case where we assumed a cooling time model. Furthermore, we daWn ”* ^ C MC ° f the cooiin S t,me model. 

again have the situation where the largest residuals comes from the two outermost points. The IPC part of the fit gives a contribution to * 2 of 26 for 26 data points which 

Therefore, because of the problems with the outermost mass points, it is not clear for our indicates a rather good fit. As in our discussion of the cooling time model, this does 


indeed represent a relatively good fit since the issues that made the HR1 data difficult to ‘han about 1.2 is good, but the model significantly overpredicts the emission at 0.65 keV 

handle are not very important in this case. This point has been discussed more fully above and underpredicts the emission at 1 keV. This is again understandable for the following 

and in section lll.(b) of TB. Finally, we move on to the FPCS part of the fit. From panel reasons. First, for the the energy range around 0.65 keV, although all the FPCS lines are 

(g) of Figure 24, we see immediately that the FPCS residuals are large and the fit to the n °‘ weU fit . the 0vm L X“ Une in particular is well fit (see panel (g) of Figure 24). It 
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resulted in a decrease in the total * 2 a* well as an increase in the abundance of oxygen. the predicted 7 keV equivalent width, we find, as in the previous cases, that an abundance 

The indicated value of * a would now seem to imply a somewhat reasonable fit (2<r) but as of more than 2 times solar is required. Additionally, the required abundance of oxygen 

before, problems with this fit as well the question of parts of the fit are improved indicate ia more than 40 times the solar value. This then makes the fit unrealistic because of the 

that no real improvement is obtained. unreasonable required abundances. 
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ft constant value &t aniftll radius. Now, in the fitting process, as before, we vary A^ali and from the ambient phase anyway, a slightly larger j> in the case of the flattened power law 

n.,. 1 , along with the other parameters of the fit but we fix o_,h for a given optimisation. will aoi result in any extra emission above the case with the power law model. Now, as for 

We then take another value for q..,, and we perform another fit to the data. This process the reason why p for the flattened power law ends up higher than in the case of the pure 

is repeated for a variety of values for P ower must look at the innermost points of the HRI where there is a significant 
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computation of the mus accretion rate but in the present context of a multiphase model 

the situation becomes considerably more complex because the flow of heat between the We end by emphasizing that our present mass accretion rate determinations represent 

ambient rpedium and the cooling blobs becomes difficult to model, improvements on the previous work that we have reviewed in the previous paragraph. The 
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assumed column density by the amount quoted hy White et al. and maJking a fit to the the SSS and the FPCS are observing basically the same gas (which is the case in our fits), 

data assuming the cooling time model. We found that although the simultaneous fit to the and if the FPCS line flux in a given energy range is well fit, then so too must the SSS 

HRI, IPC, and FPCS was indeed improved (almost all of the improvement came in the HRI spectrum in that same energy range. As it turns out, the constraint from the OVIII Lya 

part of the data), the requirement that the predicted equivalent width of the 7 keV Fe line line as observed by the FPCS implies that the emission at an energy of around 0.65 keV 

complex must also fit the measured value of about 1 keV implied unreasonably high metal in the SSS be about twice as high as is observed. This discrepancy cannot be understood 
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in the course ofithe work on the extra absorption, we were able to understand the the means to further address it here. The fact of this potential problem with the SSS 

lack of a good fit to the SSS of our models as an apparent inconsistency between the should cast considerable doubt upon conclusions that are based solely upon SSS data, for 

measurements of the SSS and the FPCS. This conclusion came from the argument that if example the ciajm of extra absorption in the M87 cluster gas. We emphasize, however, 
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FIGURE 7: FPCS rcsidu&ls. Wc plot here the residuals (the observed values m inim 

the predicted values divided hy the error) due to the FPCS lines. Line numbers on the 
abscissa correspond to observed FPCS lines and line blends. Number 1 corresponds to 
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light solid line corresponds to the the FeXVIJ line. The other lines correspond to the FIGURE 24: FPCS residuals. Panel (g) gives the FPCS residuals for the best fit 

various Fe L blends. power law multiphase model with Galactic absorption. Panel (h) gives the case where the 

column density is increased to 1.75 x 10 2, cm 3 . 

FIGURE 18: SSS spectrum. These are the predicted SSS spectra in the cooling 

time multiphase fits to the data excluding that of the FPCS. The light solid histogram FIGURE 24: p profiles. The solid line gives p for the best fit cooling time multi- 
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FIGURE 23: IPC surface brightness distribution. The dark solid line corresponds 

to the best fit power law multiphase model with Galactic absorption. The dotted line is 
for the best fit power law multiphase model with column density of 1.75 x 10 a, cm 3 . 
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